Nonuniqueness of the Potentials of Spin-Density-Functional Theory 
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It is shown that, contrary to widely held beliefs, the potentials of spin-density-functional theory 
(SDFT) are not unique functionals of the spin densities. Explicit examples of distinct sets of poten- 
tials with the same ground-state densities are constructed, and general arguments that uniqueness 
should not occur in SDFT and other generalized density-functional theories are given. As a conse- 
quence, various types of applications of SDFT require significant corrections or modifications. 
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The Hohenberg-Kohn (HK) theorem j|, which guar- 
antees that the ground-state single-particle density on its 
own is sufficient to uniquely determine all observables of 
a many-body system, is one of the most remarkable the- 
orems of quantum mechanics. It is also at the heart of 
density- functional theory (DFT), one of the most pop- 
ular many-body methods [0Jq]. In the present paper a 
fundamental problem in the extension of this theorem to 
spin-polarized situations is analysed, and consequences 
for practical applications of SDFT are pointed out. 

In the case of the original formulation of DFT, in which 
the basic variable is the particle-density n(r), the HK 
theorem can be cast in the form of two logically inde- 
pendent one-to-one maps Q . Quantum mechanics guar- 
antees that, for given particle-particle interaction U and 
particle number N, the potential in which the particles 
move determines the ground-state many-body wave func- 
tion (via solution of Schrodinger's equation), which in 
turn determines the ground-state density (by simple in- 
tegration). The essence of the original Hohenberg-Kohn 
theorem for nondegenerate ground states is that both 
of these maps are invertiblc: the ground-state density 
n(r) uniquely determines the ground-state wave function 
^(ri, . . . r^v), which in turn determines, up to an additive 
constant, the potential & 



v(r) -4=k> *(r 



• rjv) 



n(r). 
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Map two implies that every ground-state observable is a 
unique functional of the density, whereas the combined 
map allows the much stronger statement that the poten- 
tial (and, given the interaction, thus every observable) is 
a unique density functional. These abstract maps have 
found extremely important practical applications in the 
form of the Kohn-Sham (KS) formulation of DFT @, 
which is used for almost all band-structure calculations 
in solid-state physics, and a rapidly increasing number of 
electronic-structure calculations in quantum chemistry. 



Many of these applications, however, do not employ the 
original formulation of DFT, but spin-density- functional 
theory ||, in which the fundamental variables are the 
spin-resolved particle densities n CT (r). In SDFT the map 
from spin densities to wave functions is easily established, 
but that from wave functions to potentials could, in spite 
of considerable effort ^HQ] not be proven and remains 
an, albeit popular, conjecture. In the early days of SDFT 
von Barth and Hedin || already pointed out that the 
uniqueness of the spin-dependent potentials is not guar- 
anteed, and explicitly constructed two potentials which, 
when used in a one-body Hamiltonian, have common 
eigenstates. The question whether a similar construction 
is possible in the many-body case, however, remained 
open, and even in the one-body case it was objected that 
these common eigenstates were not necessarily common 
ground states ^,0. 

In this paper the general problem of uniqueness of the 
SDFT potentials is solved by (i) giving a general char- 
acterization (not limited to SDFT) of possible sources of 
nonuniqueness in DFT, (ii) explicitly constructing differ- 
ent SDFT potentials with a common ground state, and 
(iii) pointing out the physical reason behind nonunique- 
ness. 

To address point (i), let v' =v + Av and B' = B + AB 
be electrostatic and magnetic fields that are supposed to 
give the same (many-body) ground state 'F as v and B. 
Aw and AB then necessarily satisfy 
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n(r)Au(r) - m(r)AB(r)] >F = A£*. 



(2) 



A necessary condition for nonuniqueness is thus that, 
given 'F, one can find a linear combination of the density 
operators (h and m in the case of SDFT) of which ^ is 
an eigenfunction. This is automatically the case if a lin- 
ear combination exists that is a constant of motion, since 
constants of motion commute with the Hamiltonian and 
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thus have the same eigenf unctions (assumed nondegen- 
erate for convenience). If the Hamiltonian expressed in 
terms of v and B has a gap between its ground state and 
its first excited state, one can always make sure that \I/ 
remains ground state of the new Hamiltonian, contain- 
ing v' and B', by making the coefficients of the linear 
combination, Aw and AB, sufficiently small. 

From any extensive constant of motion (i.e., one which 
is a linear combination of the density operators) one can 
thus systematically construct families of potentials with 
the same ground state. As an example consider a system 
with an energy gap, for which the total magnetization 
M z = J d 3 rm 2 (r) is a constant of motion. The choice 
Av = 0, AB = Bu z — const, where u z is a unit vector 
in the z-direction, then is clearly a suitable linear combi- 
nation, as long as B is not large enough to induce level 
crossings ||. The presence of a gap is crucially impor- 
tant, since it guarantees that for sufficiently small values 
of B such level crossings do not take place, so that ^ 
remains the ground state. 

It is important to emphasize at this point that the 
exchange-correlation (xc) functionals are still guaranteed 
to be unique density functionals by the second part of the 
HK theorem. Since the Hartree potential, too, is mani- 
festly a unique density functional, the nonuniqueness can 
only arise from the external potential contribution to the 
effective potentials. Practical consequences of this are 
pointed out towards the end of this paper. 

Examples for nonuniqueness obtained from constants 
of motion are below referred to as systematic nonunique- 
ness. Apart from these, Eq. (^|) may, under suitable 
circumstances, also permit solutions not systematically 
associated with conserved quantities. To construct an ex- 
ample for such accidental nonuniqueness first recall how 
SDFT describes a spin-polarized system. After choosing 
the xc functional and specifying the corresponding spin- 
dependent Kohn-Sham potential u s . CT (r), one solves self- 
consistently the spin-dependent Kohn-Sham equation, 
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' v s,a( T ) 
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and calculates the spin densities from 



n CT (r) 



N„ 



(3) 



(4) 



where iV CT , the number of particles with each spin, is 
determined by distributing the N particles over the N 
lowest eigenfunctions <j) n<7 , such that the total energy is 
minimized. If the number of filled spin up levels is equal 
to that of filled spin down levels, the system is unpolar- 
ized (iVf = iV|), while it aquires a finite spin polarization 
if they are different. The maximum possible spin polar- 
ization is realized if all particles have the same spin (say 
spin up), so that ny = n and n± = 0. 



In the case of a fully polarized system (e.g., saturated 
ferromagnetism) , the N lowest eigenvalues of the spin up 
Kohn-Sham equation are all lower than the lowest eigen- 
value of the spin down Kohn-Sham equation, so that no 
eigenfunction of the latter is occupied. The full density 
is thus obtained from the selfconsistent solution of 



-^V 2 + Us , T (r) 



</ > nT( r ) = enT^"T( r ) 



with 
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i(t) = n T (r) =^0; T (r)0„ T (r). 
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In this case the spin up Kohn-Sham potential of SDFT 
already determines the full particle density n(r). The 
spin down potential, on the other hand, is a completely 
arbitrary function of r, as long as its lowest eigenvalue eoj 
is higher than the iV'th eigenvalue of the corresponding 
spin up equation, e_/vf (which is precisely the condition 
for the system to remain fully polarized) Q. The set of 
spin densities {rt j = n, nj, = 0} thus does not uniquely 
determine the set of potentials {w s |,u s |}, because there 
trivially is an infinite number of possible spin down po- 
tentials. Since spin densities and wave functions are still 
in one-to-one correspondence, this implies that the po- 
tentials of SDFT are not determined completely by ei- 
ther the KS Slater determinant $ or the full many-body 
ground state \&. 

In the more general case of a partially polarized sys- 
tem the above simple argument for nonuniqueness breaks 
down, but it is not hard to show that partial polarization 
does not restore uniqueness. Consider a given representa- 
tion of the ground state KS Slater determinant in terms 
of single-particle orbitals (j> n - In ordinary DFT, based 
on the density only, each of these orbitals determines the 
corresponding KS potential (up to the irrelevant additive 
constant e„ ) according to 



v s (r) 



i(j)n(r) 
<M r ) ' 
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where t is the single-particle kinetic energy operator. Eq. 
([?]) is simply the Kohn-Sham equation solved for v s (r). 
If one now tries to repeat this reasoning in the spin- 
dependent case, one must solve the Kohn-Sham equation 
of SDFT for the external potentials f s ,T(r) and v a< i(r) in 
terms of the single-particle orbitals </>„f (r) and </>„|(r), 



t(f>na(r) 



(8) 



Here we restricted us to collinear spin configurations, 
so that the magnetic field can be expressed in terms 
of the spin-dependent potentials as B z (r) — [v s ^(r) — 
u Sj |(r)]/ ji, where /i = qh/(2mc) is the Bohr magneton. 
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Again, any of the single-particle orbitals determines 
the corresponding potential up to a constant, which re- 
mains unspecified. But, unlike in the previous case one 
can choose only one constant freely by adjusting the zero 
of energy. That is, given 0„|(r), say, one can determine 
i> Sj j(r) up to an irrelevant constant, but then the zero 
of energy has been fixed, and </>n|(r) determines w s .j,(r) 
only up to a (now physically relevant) constant, rigidly 
shifting spin up levels with respect to spin down levels. 
Consequently, the KS Slater determinant does not com- 
pletely determine the spin-dependent potentials. By ap- 
pealing to the second part of the HK theorem one readily 
concludes that the same holds for the many-body wave 
function V?. 

The only caveat to the above construction is that for 
the original KS determinant to remain the ground state 
after adding the free constant to v s ,j,, this shift must 
not induce the presence of new (previously unoccupied) 
single-particle orbitals in the Slater determinant. As 
above, this is guaranteed if the KS N-particle ground 
state is separated from the first excited state by a gap, 
and the shift of w s j is sufficiently small compared to this 
gap. 

Adding a constant term, v, to v S; |(r) amounts to 
adding v to the electrostatic potential v s (r) and to 
the magnetic field B z (r). The nonuniqueness associated 
with the fact that there is only one zero of energy, but 
two constants to be determined in Eq. (||), has thus 
brought us back to a particular example of systematic 
nonuniqueness, in which the linear combination in Eq. 
(||) is given by Av — v and AB = (w//i)u z . In general we 
expect that systematic nonuniqueness is the only generic 
source of nonuniqueness, while accidental nonuniqueness 
requires certain 'pathological' features, such as complete 
spin polarization. 

Reflecting on the previous examples, what all cases of 
nonuniqueness, accidental or systematic, have in com- 
mon is that the wave functions do not change if certain 
changes are made to the potentials. This implies a de- 
gree of 'rigidity' or 'incompressibility' of the former with 
respect to the latter. The relation to incompressibility 
is, in fact, more than a mere analogy: an essential fea- 
ture of all nontrivial examples of nonuniqueness (i.e., all 
except the additive constant one can choose in the elec- 
trostatic potential) is a gap between ground and first ex- 
cited state, as a consequence of which sufficiently small 
changes to the potentials cannot bring down an excited 
state to have lower energy than the ground state. The ex- 
istence of a gap between ground and first excited state, 
on the other hand, is tantamount to incompressibility 
in the usual (quantum liquid) sense of the word. Such 
incompressibility thus constitutes a necessary (but not 
sufficient) condition for nontrivial nonuniqueness. 

In the remainder of this paper we look at several ap- 
plications of SDFT in order to determine what changes 
become necessary in view of the above results. First of 



all, in the standard Kohn-Sham formulation of SDFT the 
nonuniqueness of v s ^ a is not critical, since it arises from 
the external potentials (which are known from the begin- 
ning), whereas the xc potentials remain unique. Fields 
in which significant changes are required if the first map 
is no longer available, include (i) the calculation of ex- 
cited states, (ii) the determination of exact KS poten- 
tials from numerically exact densities, (iii) perturbative 
strategies to construct approximate density functionals, 
(iv) the optimized-effective potential method of SDFT, 
and (v) several recent generalizations of DFT. 

(i) Let us consider excited states first. A standard 
argument Q to justify the calculation of excited states 
in DFT (albeit not by direct identification of the Kohn- 
Sham eigenvalues with true excitation energies) is that 
these, too, are functionals of the ground-state density; af- 
ter all, according to the combined map w[^[ri]] this density 
determines the full potential, which, via Schrodinger's 
equation, determines all eigenstates, ground state as well 
as excited ones. Clearly, this argument is only correct in 
the original formulation of DFT, in which the combined 
map is available. The above discussion shows that it can- 
not be used to justify attempts to employ SDFT for the 
calculation of excited states [jlO| . 

(ii) Recently several methods for constructing accurate 
KS potentials, using as input highly precise numerical 
densities (obtained from Monte Carlo or configuration 
interaction calculations), have been proposed [pl]-p7[. 
Uniqueness of the resulting KS potentials is taken for 
granted in these prescriptions. Again, this holds only in 
the original formulation of DFT, but not in the (much 
more widely used) SDFT. The generalization of these 
methods to determine the KS potentials of SDFT must 
thus be reconsidered. Even the xc potentials, which are 
unique, are not completely free of this problem: in an 
exact calculation the nonuniqueness of the KS poten- 
tials must of course cancel in the final result, but this 
is not guaranteed to be the case in an actual (necessar- 
ily approximate) implementation, e.g., using finite basis 
sets. Special care must thus be taken in such calcula- 
tions to ensure that the xc potential is not contaminated 
by nonuniqueness. 

(iii) A useful method for obtaining LDA-type approx- 
imations for the xc functional is based on perturbative 
diagrammatic calculations of the xc energy of the elec- 
tron gas [^],[l8| • The results of such calculations are nec- 
essarily obtained in terms of the potentials acting on the 
unperturbed system. As long as the first map of the HK 
theorem holds, they are thus automatically also function- 
als of the densities. Obviously this ceases to be true if 
that map is no longer available. It follows from the above 
that the SDFT counterpart of this procedure works only 
because the electron gas does not have a gap between its 
ground and first excited state. Attempts to apply this 
procedure to systems with a gap (such as the supercon- 
ducting electron gas considered in Ref. [p|) thus face 



3 



unanticipated and fundamental problems in their gener- 
alization to the spin-resolved case. 

(iv) In the optimized-effective potential (OEP) method 
one minimizes an energy functional that is an explicit 
functional of the single-particle orbitals and only an im- 
plicit functional of the densities. In order to calculate the 
functional derivative with respect to the densities, one 
then employs the functional chain rule, to write |ly,K| 
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v xc [n}(r) 



SE xc [{</> n }] 



J dr' J dr" 



5n(r) 

SE xc [{</> n }] 84>n(v') 5v{v") 
S(j> n {r') 5v{r") Sn(r) 



+ c.c. 



(9) 



The first functional derivative on the right-hand side can 
be calculated explicitly, the second is easily found from 
perturbation theory, and the last is usually identified 
with an inverse response function. In SDFT the poten- 
tials are not unique, so that this inverse does not exist 
for all densities, and hence the solution to the OEP inte- 
gral equation following from Eq. (^) is not guaranteed to 
be the correct v xc . The problem here arises from using 
the KS potential as an intermediate step in performing 
the functional derivative of the orbital functional, the xc 
potential itself is of course still unique. This intermedi- 
ate step can be avoided by directly using 5(f> n /Sn on the 
right-hand side, but the resulting equation will be more 
complicated to solve than the standard OEP equation. 

(v) Finally, consider generalized density-functional 
theories. The fact that even in the oldest and most 
widely used of all generalizations of DFT, namely SDFT, 
uniqueness of the potentials is not guaranteed, casts 
serious doubts on uniqueness claims in other, more 
recent, generalizations of DFT. It can be suspected 
that in, e.g., relativistic DFT, time-dependent DFT, 
polarization-dependent DFT, DFT for superconductors, 
etc, similar problems to those discussed here in the con- 
text of SDFT are lurking. Indeed, in current-DFT J|| 
this is known to be the case, and examples of both sys- 
tematic and accidental nonuniqueness can be found p2[ . 

In summary, it has been shown, both by explicit ex- 
amples and by general considerations, that, unlike in or- 
dinary DFT, in SDFT and other generalized DFTs the 
effective and the external potentials are not uniquely de- 
termined by the spin densities alone. Nonuniqueness can 
arise accidentally, via special features of the ground state, 
and systematically, via extensive constants of motion. In 
both cases it depends on the presence of a gap in the 
spectrum of the system, and on the resulting incompress- 
ibility. As a consequence many previous applications of 
SDFT and other generalized DFTs in condensed matter 
physics and quantum chemistry must be critically reex- 
amined. 
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